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ABSTRACT 



Context. Gamma-ray burst (GRB) afterglows are excellent and sensitive probes of gas and dust in star-forming galaxies at all epochs. 
It has been posited that dust in the early Universe must be different from dust at lower redshifts. To date two reports in the literature 
directly support this contention, one of which is based on the spectral shape of the afterglow spectrum of GRB 050904 at z = 6.295. 
Aims. Here we reinvestigate the afterglow of GRB 050904 to understand cosmic dust at high redshift. We address the claimed evidence 
for unusual (supernova-origin) dust in its host galaxy by simultaneously examining the X-ray and optical/near-infrared spectrophoto- 
metric data of the afterglow. 

Methods. We derive the intrinsic spectral energy distribution (SED) of the afterglow at three different epochs, 0.47, 1.25 and 3.4 days 
after the burst, by re-reducing the Swift X-ray data, the 1.25 days FORS2 z-Gunn photometric data, the spectroscopic and z'-band 
photometric data at ~ 3 days from the Subaru telescope, as well as the critical UKIRT Z-band photometry at 0.47 days, upon which 
the claim of dust detection largely relies. 

Results. We find no evidence of dust extinction in the SED at any time. We compute flux densities at /t rcst = 1250 A directly from the 
observed counts at all epochs. In the earliest epoch, 0.47 days, where the claim of dust is strongest, the Z-band suppression is found to 
be smaller (0.3 ± 0.2 mag) than previously reported and statistically insignificant (< 1.5cr). Furthermore we find that the photometry 
of this band is unstable and difficult to calibrate. 

Conclusions. From the afterglow SED we demonstrate that there is no evidence for dust extinction in GRB 050904 - the SED at all 
times can be reproduced without dust, and at 1 .25 days in particular, significant extinction can be excluded, with A(3000 A) < 0.27 mag 
at 95% confidence using the supernova-type extinction curve. We conclude that there is no evidence of any extinction in the afterglow 
of GRB 050904 and that the presence of supernova-origin dust in the host of GRB 050904 must be viewed skeptically. 

Key words. Galaxies: high-redshift - ISM: dust, extinction - Gamma rays: bursts 
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1. Introduction 

Gamma-ray burs ts (GRBs) can be detected up to the ons et of 
reionization (e.g. iTanvir et al]|2009t ISalvaterra et al.ll2009h due 
to their brightness in th e first few hours after the explosion 
dLamb & Reichartll2000l) . GRBs are transient sources followed 
by long lasting afterglows, emitting energy intensely across the 
full range of the electromagnetic spectrum. GRB afterglows are 
effective and informative probes for the detailed study of cosmic 
dust at high redshifts due to their simple power-law spectra, high 
luminosities and locations in star-forming environments. 

Interstellar dust has a crucial significance in the appear- 
ance and evolution of star formation in the early Universe. 
It is still under debate whether interstellar dust properties 
have evolved as a function of redshift. At high redshift {z 
>5-6) it has been suggested that dust might originate in 
sources other than the evolved asymptotic giant branch stars 
that are the dominant source of dust produ ction in the lo- 
cal Universe dGehrzl Il989t iDwek et all I2007I) . Previous stud- 
ies reported that dust in cosmological objects at z > 6 is 
predominantly produced in the ejecta of core-collapse super- 
novae (SNe), rath er than the evolved stars which are missing on 
short timescales dTodini & Ferraral l200lb iMorgan & Edmundsl 



20031: iNozawa et al.l|2003l: IDwek et alj | 2007t lMarcherikfll[ 200e 
Hirashita et al.ll2005l) . Recently, however. IValiante et alJ (2009 



calculated that the most massive asymptotic giant branch (AGB) 
stars could produce dust on time scales sh ort enough to dominate 
dust production by z ~ 6. Observationallv. lMaiolino et al.l d2004h 
found an unusual extinction curve in the most distant known 
broad absorption line quasar (BAL QSO) at redshift z = 6.2, 
consistent with what could be expected from dust produced in 
core-collapse SNe. 

The progenitors of long -duration GRBs are known to 
be short-lived massive stars dGalama et al. 1998t iHjorth et al.l 
l200lRIStane"k et al.l2003l:lMalesani et alJ2004h . GRB 050904 at 
z = 6.295 was a long duration GRB. It was extremely luminous 
and is the third most distant known GRB to date. GRB 050904 
was detected by Swift on 2005 September 4 at f = 01:51:44 UT 
dCummings et al.l 120051) . Substantial multi-wavelength follow- 
up was carried out simultaneously for GRB 050904 with sev- 
eral ground based facilities. Previous an alysis of the rest-frame 
UV afterglow found no evidence of du st ( Tagliaferri et alJl2005b 
Haislipet al l 120061: iKann et al.l 120071: iGou et all 12007^ Later 



Stratta et al.l d2007 ) re-examined the afterglow SED at different 
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epochs and claimed the detection of dust absorptiontion with an 
extinction curve consistent with that used to explain the spec- 
trum of the highest redshift BAL QSO, but inconsistent with the 
dust typical of the local Universe. 
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The claim of detection of SN-origin dust in GRB 050904 is 
of fundamental importance to the question of the origin of dust 
in the early Universe, a very vexed proble m for high redshift 
sub-m m galaxies (see, e.g. the discussion in iMichalowski et all 
It was the first only direct observational evidence of dust 
from SNe in a high redshift star-forming environment. In this 
paper we review the relevant data to test whether dust is required 
by these observations and if so, what kind of dust is needed. The 
outline of the paper is as follows: In §2 we describe the detailed 
multi-band spectral analysis of the afterglow of GRB 050904 at 
different epochs. In §3 we present results from the SED of the 
afterglow. In §4 we discuss possible scenarios. In §5 we provide 
our conclusions. 

2. Multi-wavelength observations of the afterglow 

2. 1 . X-ray analysis 

Swift's, X-ray Telescope (XRT) localized the afterglow of 
GRB 050904. The XRT data (in the 0.3-10.0 keV energy range) 
were extracted and reduced using the HEAsoft software (version 
6.4). We computed the X-ray spectra at three epochs, specifically 
0.47, 1.25 and 3.4 days post-burst, chosen as the epochs with the 
best spectroscopic and photometric optical/near-infrared cover- 
age. X-ray spectra at three epochs were created in a standard way 
using the most recent calibration files. 

For our analysis, it is important to obtain an accurate es- 
timate of the absolute flux for these X-ray spectra. The X- 
ray lightcurve i s extremely variable, exhibiting long lasting 
flaring activity (|Watson et al.l l2006bt ICusumano et al.l 120061 : 
iGendre et al.l 120061) . The flares suggest two separate compo- 
nents, which may be due to a num ber of causes, possib ly activ- 
ity of the GRB central engine (e.g. lBurrows et al.ll2007h . At late 
times the X-ray count rate is very low (see Fig.|2|, therefore, we 
aim to get an accurate estimation of the flux which includes the 
uncertainty due to the flares. We fitted the aft erglow lightcurve 
by as suming a smoothly broken power-law dBeuermann et al] 
1 19991) to get an approximate overall X-ray lightcurve. The fit 
to the afterglow lightcurve results in at most a 30 - 40% scatter 
around this fit at all epochs. We then normalized the X-ray spec- 
tra to the lightcurve fit at the relevant SED epoch. The procedure 
results in X-ray spectra with the best estimate of the slope and 
flux at the relevant SED epoch, and with an additional overall 30 
- 40% uncertainty on their absolute flux levels. 

2.2. Near-IR and optical imaging 

Several telescopes obtained photometric obser vations of the af- 
terglow in the optical and near-infrared ban ds (ITagliaferri et al.l 
120051: lHaislip et al.ll2006fc iPrice et al.ll2006l) . For a comprehen- 
sive investigation of the SED, we g athered optic a l and near- 
infrared photometry at three epochs. IStratta et alJ (120071) sug- 
gested unusual dust particularly on the basis of the Z-band ob- 
servation at ~ 0.5 days, taken by the UK Infrared Telescope 
(UKIRT) equipped with WFCAM. Therefore, we re-reduced 
these Z-band data obtained at fo + 0.47 days. The object has 
low signal to noise. Our best estimate of the magnitude comes 
from point-spread function (PSF)-matched photometry, carried 
out using the DAOphot package within IRAF. For consis- 
tency we also carried out aperture photometry (with Gaia and 
IRAF/PHOT), and found the afterglow Z-band magnitude to be 
strongly sensitive to the chosen sky extraction annulus, being in 
some cases brighter by ~ half a magnitude, although with large 
errors. This uncertainty is worth noting. 



At 1.25 days z-band photometry was obtained with the 8.2m 
ESO Very Large Telescope (VLT) by using the FOcal Reducer 
and low dispersion Spectrograph 2 (FORS2) z-Gunn filter. The 
object is detected with high signal to noise and using aperture 
photometry we could s ignificantly reduce the error reported by 
ITagliaferri etall (120051) . which was likely mostly due to calibra- 
tion issues. 

At 3.3 days after the burst, z'-band phot ometry was car- 
ried out with the 8.2m Subaru telescope dive et al.l 120041) 
using the Faint Objec t Camera and Spectrograph (FOCAS; 
Kashikawa et a I] 120021) . We re-analysed the photometric data 
using aperture photometry. We compared our FOCAS z' re- 
sult with the Gemini South G MOS z' value at 3.16 days ob- 
tained by Ha islip et al.l d2006h . Correcting for the time differ- 
ence between the two observations using a broken power-law 
with ai = 0.72 + g;^, g 2 = 2.4 + 0.4, f b = 2.6 + 1.0 days 
dTagliaferri et al.l 120051) . we found these photometric measure- 
ments to be consistent within the statistical uncertainties. 

Other photome t ric data were taken from the literature 
dHaislip et al.ll2006h ITagliaferri et al.ll2005l) when available at a 
time close to our nominal SED epochs (/-band at 0.49 days). 
All H and K s data as well as J at 1. 25 and 3.4 da y s wer e de- 
rived from the best-fit lightcurves of lHaislip et al.l d2006l) and 
ITagliaferri et al.1 (120051) . At 1.25 days we used the F-band pho- 
tomet ry from the lightcurve computed in this band (lHaislip et al.l 
12006b . 

We corrected the observed magnitudes for extinc tion in 
the Milky Way (MW) using the ISchlegel etall dl998l) maps 
with R v = 3.1 and E(B - V) = 0.081 mag along the line 
of sight of the burst. Potential systematic uncertainties in the 
Galactic extinction correction have no significant effe ct on our 
res ults. The independent a nalysis of lDutra et al.l d2003l) confirms 
the ISchlegel et all d!998l) maps for low E(B - V). For the z- 
bands considered here, even assuming a 30% uncertainty in the 
Galactic value would correspond to an extinction uncertainty of 
0.03 mag which is smaller than the statistical uncertainties we 
find for the extinction in the host galaxy (see Table[T]i. Effects in 
the J, H, and K s bands will be smaller still. An under-estimate 
of the Galactic extinction would lead to a smaller host galaxy 
extinction than we find in the current analysis. 

The z-Gunn, FOCAS z', GMOS z' and Sloan Digital Sky 
Survey (SDSSB iFukugita et all 1 19961) z filters have almost the 
same profile across the whole band and extend much redder than 
the UKIRT Z filter. It should be noted that since the Lya absorp- 
tion occurs in these bands, the filter wavelength coverage affects 
the observed magnitude significantly (see § 3.2). Unless explic- 
itly mentioned, in the rest of the paper we use the term "z-band" 
to denote all of the three z filters, i.e. UKIRT Z, FORS2 z-Gunn 
and FOCAS z'. 

2.3. Grism spectroscopy 

The optical spectrum of the afterglow was obtained with FOCAS 
at the Subaru telescope and the spectroscopic dat a were re- 
trieve d from the SMOKA science archive facilitjQ (Babaet al. 
120021) . The afterglow was observed on 7 September 2005 for 
4.0 hr. The exposure mid-time was 12:05 UT, corresponding to 
3.4 d ays after the burst trigger dKawai et al.l 120061 iTotani etall 
120061) . The spectra were flux calibrated using the spectrophoto- 
metry; standard star BD+28D4211 obtained on the same night. 
Individual spectra were combined following standard data re- 

1 http://www.sdss.org/ 

2 |htt p : //smo ka . nao . ac . jp/1 
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duction techniques using IRAF. The spectrum exhibits no flux 
below ~ 8900 A, consistent with a break due to Lya absorp- 
tion at redshift z ~ 6.3 and the Lya forest. The spectrum shows 
a flat continuum at the red wavelength end, revealing a series 
of metal absorption lines arising from different atomic species 
at z = 6.295, and an intervening C iv system at z = 4.84 
dKawai et al.ll2006h . The observed spectrum was corre cted for 
Galactic extinction by assuming the lCardelli et al.l d 1989b extinc- 
tion curve and as explained in § 2.2 above. 

We implemented Voigt profile fitting to the 3.4 day 
Subaru spectrum using th e FITLYMAN package in MIDAS 
dFontana & Ballesterl [19951) . We measure a hydrogen column 
density of logA^Hi (cm' 2 ) _ 2 1 .62 + 0.02, consistent with the 
value reported by iTotani et"aT1 (120061) . It should be noted that 
z'-band photometry and spectroscopy of the afterglow were ob- 
tained with FOCAS at 3.3 and 3.4 days, respectively. 

3. SED analysis 

IStratta et all d2007l) studied the optical-UV rest-frame SED of 
the afterglow of GRB 050904 at 0.5, 1 and 3 day epochs and 
found a deficit in the z-band at 0.5 and 1 days, and (less signifi- 
cantly) at 3 days, compared to the JHK power-law extrapolation, 
claiming that dust reddening could explain the flux deficit. This 
required a SN-type extinction curve. 

3.1. Afterglow compound SED 

Knowledge of the SED can address the z-band flux suppression 
issue, therefore, we computed the near-infrared to X-ray SED 
of GRB 050904 at three epochs, i.e. 0.47, 1.25 and 3.4 days. To 
facilitate comparison of the z-band flux, the SED at all epochs 
was normalized to the //-band flux, using the smoothly broken 
power-law presented bv lTagliaferri et all d2005h . The normalized 
near-infrared photometry is generally consistent, but the X-ray 
spectra are much brighter at 0.47 and 1 .25 days due to the in- 
tense afterglow flaring activity at these times. The consistency of 
the X-ray flux with the NIR SED extrapolation suggests that the 
X-ray afterglow at 3.4 days was relatively unaffected by flares. 
The composite SED of the afterglow of GRB 050904 at three 
different epochs is shown in Fig. [2] 

3.2. Comparing the z-band filter responses 

Since the z-band photometry is strongly affected by the Lya ab- 
sorption, we performed spectro-photometric analysis by utiliz- 
ing the total effective filter transmission functions including de- 
tector responses (Fig. [TJ. We use the following method that al- 
lows for a clean comparison of the different z-band magnitudes 
of the afterglow, taken at 0.47, 1.25, and 3.4 days after the burst, 
with the filters UKIRT Z, VLT z-Gunn, and Subaru z', respec- 
tively. The method essentially constructs the SEDs of stars in the 
field and uses these to make a direct comparison of the afterglow 
magnitudes at each epoch. 

First, in each afterglow image we select several non- 
saturated reference sta rs with known SPSS a nd Two Micron All 
Sky Survey (2MASS3; Skrutskie et alJl2006h magnitudes. Using 
the 2MASS /-band, and the SDSS z and i, we construct a rudi- 
mentary SED for each reference star, where we convert the mag- 
nitudes to flux densities (in erg s crrT 2 A~') at the central 
wavelength of the SDSS and 2MASS filters. We connect these 

3 http : //www . ipac . caltech.edu/2mass/ 
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Fig. 1. Filter transmission curves of SDSS z, FOCAS z', FORS2 
z-Gunn and UKIRT Z. The grey curve corresponds to the 
median-filtered optical spectrum at 3.4 days (arbitrarily normal- 
ized). The vertical triple black dot-dashed line represents the 
wavelength (1250 A) where we computed the flux density (see 
§3.2). 

flux densities with a linear interpolation, and integrate the ref- 
erence star SED convolved with the filter response curve rele- 
vant to that image, retrieving the band-integrated flux in erg s _1 
cm -2 . Second, from the image we measure the counts for the 
reference stars using aperture or PSF photometry, determining 
the conversion factor between counts and flux. Using this fac- 
tor, we eventually compute the (band-integrated) afterglow flux 
from its measured counts. We used several comparison stars to 
evaluate the accuracy of the procedure. At 0.47, 1.25, and 3.3 
days, we find a scatter of 0.02, 0.04, and 0.02 mag using 8, 10, 
and 5 reference stars, respectively. The small scatter confirms the 
robustness of our method. 

In the deep FORS2 and FOCAS z-band images, the bright- 
est stars are heavily saturated, and suitable reference stars are 
lacking since the fainter stars have large uncertainties in the 
SDSS and 2MASS catalogues. Therefore, we calibrated a set 
of faint stars using the UKIRT J and z-band images, based on 
the 2MASS and SDSS catalogs. For the z band, due to the dif- 
ference in the UKIRT and SDSS filters, appropriate color terms 
were taken into account, achieving a photometric accuracy of 
« 0.02 mag. Given the higher sensitivity of the SDSS in the i 
band, suitable calibrators for the VLT and Subaru images were 
available directly from the SDSS catalog. Note that our calibra- 
tion is entirely based on the SDSS and 2MASS catalogs, there- 
fore, our analysis is not dependent on the sky conditions when 
the data have been acquired. 

The third and final step is to convolve the relative afterglow 
spectral shape (as measured from the Subaru spectrum that was 
obtained at 3.4 days) with the three different z-band filter re- 
sponse curves, where the spectrum is rescaled in absolute terms 
to recover the band-integrated flux (in erg s _1 cm~ 2 ) determined 
for each epoch (see above). We note that this method does not 
rely on the absolute flux calibration of the Subaru spectrum; 
it merely uses the photometry to rescale it, and therefore the 
errors only include the errors in the aperture/PSF photometry, 
the error from the conversion factor, and the Subaru noise error 
when convolving it with the filter response curves. After rescal- 
ing of the spectrum, the afterglow flux density (in erg s cm -2 
A~') can be compared at any pivot wavelength, after the //-band 
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Fig. 2. Near-infrared to X-ray spectral energy distribution of 
the afterglow of GRB 050904 at 0.47 (red triangles), 1.25 (blue 
squares) and 3.4 days (black circles) after the burst. The SED 
at 0.47 and 1.25 days is scaled to the //-band at 3.4 days. The 
solid grey curve represents the median-filtered optical spectrum 
at to + 3.4 days. The black dashed line corresponds to a power- 
law fit to the near-infrared to the X-ray data at 3.4 days with a 
spectral index /3 = 1.25 + 0.05. The grey shaded area represents 
the lcr uncertainty on the power-law. 



normalization. At all epochs we computed the flux density at 
/Irest = 1250 A, which was selected since it is close to the peak 
of all the involved filter transmission curves (see Fig. [T}, and is 
separated from the metal absorption lines visible in the spectrum 
dKawai et al . 200g). 

The single but important assumption in this method is that 
the spectral shape of the afterglow is not changing from 0.47 to 
3.4 days. This is in some sense the null hypothesis that we are 
trying to test: dust destruction would produce a change in the 
relative spectral shape, and would therefore produce a change 
in the z-band brightness relative to the //-band normalization. 
Other effects, such as a variable H I column, or a change in the 
spectral slope due to e.g. the cooling frequency crossing the z- 
band, could also in principle cause such a change. However, if 
the resulting afterglow z-band brightness (normalized to the H- 
band) between 0.47 and 3.4 days does not show a significant 
change, then this would provide strong support for the null hy- 
pothesis that the spectral shape is not changing. 

Following the above procedure, we find the afterglow to have 
a flux density at 1250 A of 14.7 + 1 .32, 9.41 +0.24 and 2. 1 3 +0.22 
//Jy at 0.47, 1.25 and 3.3 days, respectively. We find that the nor- 
malized 0.47 day UKIRT Z-band brightness is 0.27 + 0.18 mag 
fainter than the FOCAS z'-band brightness at 3.3 days. At 1.25 
days, the normalized FORS2 z-Gunn brightness is brighter by 
0.17 + 0.15 mag compared to the FOCAS z'-band brightness at 
3.3 days. The uncertainties here also include the uncertainties 
in the normalization, i.e. the errors in the //-band photometry 
(0.06 mag). Therefore, there is no evidence for variability of the 
spectral shape around the z band. In particular, after taking into 
account the appropriate filter shapes and color effects, there is no 
significant deficiency of flux in the z-band flux at 0.47 days c om- 
pared to later epochs dHaislip et alJl2006tlStratta et alJl2007l) . 



Days 


Model 


P 


A(3000A) 
(mag) 


(mag) 


0.47 


PL 
PL+SN 
PL+SMC 


1.28 ±0.11 

1.22 ±0.24 

1.23 ±0.08 


0.3 ± 0.22 
0.1 ±0.07 


0.05 ± 0.04 


1.25 


PL 
PL+SN 
PL+SMC 


1.24 ± 0.09 
1.27 ±0.2 
1.17 ±0.51 


0.05 ±0.11 
0.01 ± 0.04 


0.01 ± 0.02 


3.4 


PL 
PL+SN 
PL+SMC 


1.25 ±0.05 

1.23 ±0.21 

1.24 ±0.07 


0.22 ± 0.24 
0.056 ± 0.059 


0.042 ± 0.044 



" The SN-origin extinction curve has been only computed in the 
range /l lcst = 1000 - 4000 A, hence it is not possible to provide Ay. 



4. Discussion 

At 0.47 days post-burst, we find a flux deficit in the UKIRT 
Z-band compared to the 3.3 days Subaru photometry that is 
only significant at < 1.5<x level. This low significance result, 
combined with the difficulty in determining the Z-band magni- 
tude at 0.47 days alluded to in §2.2, suggests that a change in 
the spectrum between 0.47 and 3.4 days does not have strong 
observational support. If the effect were real then such a flux 
de ficit could be expla ined by: (i) dust extinction as suggested 
bv lStratta et al.l (l2007h with a SN-origin e xtinct ion curve, or (ii) 
gas absorption. Previously lHaislip et alJ d2006l) also suggested 
that absorption due to molecular hydrogen could give rise to the 
Z-band flux deficit at 0.47 days. 



4.1. Dust in the GRB 050904 host galaxy 

The claim of SN-type dust in GRB 050904 is important be- 
cause of the possib ility of obse r ving t he evolution of cosmic 
dust at high redshift. lStratta et al.ld2007l) suggested SN-type dust 
extinction in the host galaxy of GRB 050904 with an extinc- 
tion c urve inferred for a BAL QSO at z = 6.2 dMaiolino et al.l 
l200l . The unusual extinction curve is rather flat at longer wave- 
lengths and steeply rises at A < 1700 A. The best-fit estimates of 
IStratta et al.1 d2007l) of the extinction at 3000 A in the rest-frame, 
A(3000 A), were 0.89 + 0.16, 1.33 + 0.29, and 0.46 + 0.28 mag 
at 0.5, 1, and 3 days, respectively. 

It is clear from our broad-band SED at 3.4 days (see Fig. |2| 
that the extrapolation of the near-infrared power-law is consis- 
tent with a single power-law to the X-ray spectrum, i.e. consis- 
tent with both the slope and flux level of the X-ray spectrum at 
that time. We can also clearly see that there is no evidence in the 
flux-calibrated optical/near-infrared spectrum at 3.4 days for any 
extinction - the continuum just redward of the Lyo- absorption is 
consistent with the single JHK power-law. Both facts mean that 
there is no evidence for dust extinction at 3.4 days. We fitted 
a dust-attenuated power-law using a d ust model for the Small 
Magellanic Cloud (SMC, Rv = 2.93: |Pelll992l) and the SN- 
origin extinction model of iMaiolino et al.l d2004l) to the 0.47, 
1.25 and 3.4 day zYJHK s data (from the z-band, we compute 
the flux density at /l rest = 1250 A). The best fit parameters are 
reported in TableQ] With our r evised z-band photo metry, extinc- 
tion at the level suggested by IStratta et al.l d2007l) can be ruled 
out at all three epochs (see Fig. [3}. In no case the computed ab- 
sorption is significant at more than 1.5cr level. 
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Extinction-correct ing the 1.25 day SED at the level fitted 
bv lStratta et al.l d2007l) makes its extrapolation overshoot the X- 
ray spectrum, hinting that ~ 1 mag of extinction at 3000 A is 
not required. More importantly, the 7-band photometry, with 
a central wavelengt h of 1400 A in the rest-frame, at 1.25 days 
(lHaislip et al.1120 06). is consistent with the near-infrared power- 
law extrapolation. Such consistency would not be expected in 
the lStratta et al.l (120071) dust hypo thesis since A(1400 A) is about 
1.75 times the A(3000A) in the iMaiolino et alj d2004l) model, 
and the T-band photometry should therefore lie a factor of 2 be- 
low a power-law extrapolation, while it does not (Fig. [3), though 
its error is large. As it can be seen in the middle panel of Fig.[3j 
the SED at 1 day follows a simple power-law and provides strong 
constraints on dust absorption. Again, it seems likely that not 
only is there no evidence for SN-type extinction in GRB 050904 
after 1.25 days, but that there is no evidence for any dust extinc- 
tion at all at ~ 1 day or later. 

There are also strong arguments against a SN-origin 
dust interpretation at 0.47 days. While dust reddening has 
been uneq uivoc ally observed in lower redshift GRBs (e.g. 
Kann et alJ l2006t lFvnbo et al.ll2QQ9fc GRB 05040 1 : IWatsonetalJ 



Rest wavelength (A) 
2000 



2006a!; GRB 991 216: IVreeswiik et alj [2006t GRB 050408: 



Foley et al.l 120061: Ide Ugarte Postigo et alJ l2007t GRB 070 802: 
Elfasdottir et alj|2009l: GRB 080607: iProchaska et alj|2009l) . so 



far SN-origin dust has never been seen before in any GRB 
host. Moreover there is no compelling evidence of dust ex- 
tinction in any GRB beyond z — 5. A possible exception is 
GRB 07 1025 (whi ch has a photometric redshift 4.4 < z < 5.2 
iPerlev et"aTI l2009h . which shows indications of a significant 
dust column. Notable are the two bursts at higher redshift tha n 
GRB 050904, i.e. GRB 0809 13 at z = 6.7 (iGreiner et al.ll20Q9h 
and GRB 090423 at z = 8.2 (iTanvir et alJl2009t ISalvaterra et alJ 
l2009h . neither of which show any sign of extinction. Second, 
given that dust can be excluded at t > 1 day, having non-zero ab- 
sorption at t — 0.47 days would require time-varying dust extinc- 
tion, which has never been observed in any burst. If due to dust 
destruction, we would expect reddening variations to be associ- 
ated with intense episodes of emission, while there is no optical 
flaring or any significant feature in the restframe-UV lightcurve 
in th is period that could be responsible for such dust destruction 
(see lHaislip et alJ 120061: iTagliaferri etaD l2005h . and most dust 
destruction scenarios sublimate dust on timescales of only a few 
minu t es after the burst at most (iPerna et al ] l2003HFruchter et al. 
l200ll) . IStrattaet all (120071) suggested that varying extinction may 
also indicate that the emitting region had become larger than the 
obscuring cloud. While this cannot be excluded, such a geometry 
requires some tuning of the cloud and fireball parameters. The 
claim of dust in the host galaxy of GRB 050904, with an unusual 
extinction curve, relying principally on a smaller (0.3 mag) and 
< 2cr flux deficit in a photometric observation, is not the most 
likely explanation. The most likely hypothesis is simply system- 
atic uncertainties related to the Z-band calibration. 

It is worth noting however that time-variable dust with an un- 
usual extinction curve is not even the simplest explanation even 
if the original analysis had been reliable. Given that we know 
from the optical spectrum that a large quantity of gas is present 
in the system, a variability in the gas column density at early 
times is a less tortured hypothesis. 



4.2. X-ray absorption 

The X-ray spectral analysis suggests a high metal column 
density in the afterglow of GRB 050904 dWatson et al.ll2006bt 
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Fig. 3. Near-infrared spectral energy distribution of the afterglow 
of GRB 050904 at 0.47 {top panel), 1.25 {middle panel) and 3.4 
days {bottom panel) after the burst. The observed data are cor- 
rected for Galactic extinction (§ 2.2). The corresponding bands 
are identified in the bottom panel. The solid, dashed, and dotted 
lines represent the best fit with a power-law, a power-law with 
SN dust, and a power-law with SMC dust, respectively. At 1.25 
days, the three lines almost overlap. 



ICampana et al.l 120071) . Time-resolved X-ray spectroscopy re- 
veals that the column d ensity of metals within the first few 
hours is highly variable (ICampana et al.ll2007t ICusumano et al.l 
120061: iGendre et ai1l2006l) . Due to the rapid changes in the X-ray 
spectrum this apparently variable column may be an artifact of 
the changing intrinsic spectrum resulting in a downturn at sof t 
energies that disappears at later times (see iButler et al.l [2006b . 
However, even if the change in the soft X-rays is really due to a 
variable column density, i.e. due to increasing ionization of the 
metals, this effect occurs at early times (< 10 3 s) and cannot sup- 
port the idea of dust destruction after 0.5 days. Indeed, a varying 
metal column density at < 1000 s argues against dust destruction 
at 0.5 days. If the varying metal column density is a real effect, 
destruction of any dust associated with the high metal column 
should have been completed long before 0.5 days. As a more 
general point, the optical and X-ray fluxes are at least one to two 
orders of magnitude lower after 0.5 days than before 1000 s. It is 
difficult to construct a scenario in which significant dust destruc- 
tion occurs in the interval 0.5-3 days that did not occur before in 
the absence of a huge flare in the UV-X-ray, something which is 
not observed. 

4.3. Gas-to-dust ratio 

GRB s typically occur in host galaxies with h i gh gas to dust ratios 



(e.g. iJensen et al .1 12001 



I2003at IStratta et al.l 1200 



Galama& Wiiersl l200lt iHiorth et all 
lElfasdottir et al.N2009l) The Hi col- 
umn density of the host of GRB 050904 is very large while 
Ay is small. Using our limit on (SMC-type) dust at 1.25 days, 
Ay < 0.05 mag at 95% confidence, leads to a high gas-to-dust 
ratio N{Ui)/A v > 8.3 x 10 22 cm" 2 mag _1 . The Galactic relation 
between H 1 column density and dust reddening is j V(H i)/Ay = 
4.93/^yX 10 21 ctrr 2 mag" 1 dDiplas & Savage! 19941). Correcting 



for metallicity at 3.4 days (Z = O.IZ^: iKawai et al.ll2006h . this 
implies an N{Hi)/Ay ratio l imit 5 times the Ga lactic value. A 
comparison with the SMC dGordon et al.l 1200*31) . which has a 
similar metallicity to the environment of GRB 050904, yields a 
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gas-to-dust ratio which is also more than 5-10 times larger in the 
host of GRB 050904. 



4.4. The origin of dust in the early Universe 

In the local Universe, the major sources of interstellar dust are 
AGB stars, the lower mas s ranges of which require > 1 Gyr 
to evolve to produce dust dDwek et"ai]|2007l) . It has been sug- 
gested that for sources with large dust masses such as sub-mm 
galaxies, due to the short time available at z > 5, an alterna- 
tive source of dust is required and that c ore-collapse SNe could 

ni & Ferrarall200l|: 
iDwek et all 



dominate dust forma t ion at t hese times dTod 



Morgan & Edmunds 20031: Nozawa etal.1 



2003; 



20071: lMarchenkoll2006t iHirashita et al.ll2005l) . However, more 



complete theoretical models including dust destruction by 
supernova shock or grain growth/destruction in the inter- 
stellar medium obtain yie lds that are < 0.01 M per SN 
dBianchi & Schneider! |2007|). consistent with a lmost all obser- 
vations of nearby SNe d Wooden et aH 1993 ; lElmhamdi et al.l 
12003b iMeikle et all 120061: iBlair et all 120071: ISakon et al.ll2009l) . 
This is too little to pr oduce the quantities of dust observed at 
high redshift. Recently lValiante et al.l d2009i) argued that on short 
timescales massive AGB stars could form much of the dust, de- 
pending on the assumed initial stellar metallicity and star for- 
mation hist ory. The galaxy-SED mode lling of sub-mm-selected 
galaxies of iMichalowski et al.l d2010h suggests dust-formation 
timescales of order tens of millions of years in a few cases at 
z > 4, which would clearly preclude even high-mass AGB dust- 
formation. While intriguing, these cases may be affected by ac- 
tive galactic nuclei (AGN) contamination and must be treated 
cautiously. 

Observationally, after our analysis here of the afterglow of 
GRB 050904, the detection of a peculiar extinction c urve in a 
BAL QSO spectrum at z = 6.2 dMaiolino et al.ll2004 remains 
the only direct evidence for domina nt SN-origin dust in the early 



Universe (but see recent wo rk by Perlev et al. ( 20091) ). While 



the observational analysis of iMaiolino et al.1 ( 20041) is carefully 



done, the relatively narrow wavelength coverage, the presence of 
strong, broad absorption and emission lines that dominate over 
the continuum at the blue end of the spectrum, and the use of 
composite QSO spectra, leave the result awaiting further con- 
firmation. Furthermore, it is diffic ult to exclude tha t the dust is 
affected by the central AGN itself dPerna et al.ll2003l) . so that the 
extinction curve may not tell us a lot about the origin of that dust. 



5. Conclusions 

In this work we reinvestigated the afterglow of GRB 050904 at 
0.47, 1.25 and 3.4 day epochs to understand stellar environments 
and interstellar dust at high redshift. We find that the afterglow 
SED can be reproduced at all epochs without any dust extinc- 
tion. The previous findi ng of dust extinctio n requiring a SN- 
type extinction curve by IStratta et al.l (120071) relies mostly on a 
Z-band photometric point at 0.47 days which we find has cali- 
bration difficulties and with our new accurate analysis technique 
we find the flux deficit to be both smaller and less significant 
than reported by previous studies. We can reasonably exclude 
the presence of substantial quantities of any type of dust in this 
GRB host galaxy at all epochs. We therefore conclude that there 
is no significant evidence of dust extinction in the afterglow of 
GRB050904. 
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No te added post-submission: A recent paper bv lPerlev et al.l 
(120091) reports significa nt SN-origin dus t extinction in 
GRB 071025 at z ~ 5 dPerlev et al.l 120091) . We note that 
iPerlev etaO (120091) also independently attempted to model 
the dust profile of GRB 050904 and found that the data are 
consistent with no extinction at all. 
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